% 创建时间：2025-06-22 12:12:59
% 作者：兪馳
clc;
clear;
A=sym([-3,-1,-1; ...
    0,-3,-1; ...
    1,2,0]);
syms k,
[V,J]=jordan(sym(A))

disp('----------------------------------')

A0=-2*eye(3);
H=J-A0;
J1=A0^k+k*A0^(k-1)*H+k*(k-1)/2*A0^(k-2)*H^2,
F=simplify(V*J1*inv(V))

latex(F)